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Abstract 

In this paper we investigate the vacuum polarization effects associated with a massive 
fermionic field due to the non-trivial topology of the global monopole spacetime and boundary 
conditions imposed on this field. Specifically we investigate the vacuum expectation values 
of the energy-momentum tensor and fermionic condensate admitting that the field obeys the 
MIT bag boundary condition on two concentric spherical shells. In order to develop this 
analysis, we use the generalized Abel-Plana summation, which allows to extract from the 
vacuum expectation values the contribution coming from a single sphere geometry and to 
present the second sphere induced part in terms of exponentially convergent integrals. In the 
limit of strong gravitational field corresponding to small values of the parameter describing 
the solid angle deficit in global monopole geometry, the interference part in the expectation 
values are exponentially suppressed. The vacuum forces acting on spheres are presented as 
the sum of self-action and interaction terms. Due to the surface divergences, the first one is 
divergent and needs additional renormalization, while the second one is finite for all non-zero 
distances between the spheres. By making use of zeta function renormalization technique, 
the total Casimir energy is evaluated in the region between two spheres. It is shown that 
the interaction part of the vacuum energy is negative and the interaction forces between the 
spheres are attractive. Asymptotic expressions are derived in various limiting cases. As a 
special case we discuss the fermionic vacuum densities for two spherical shells on background 
of the Minkowski spacetime. 



PACS number(s): 03.70.-Fk, 04.62.-hv, 12.39.Ba 



1 Introduction 

It is well known that different types of topological objects may have been formed in the early 
universe after Planck time by the vacuum phase transition i2| . Depending on the topology of 
the vacuum manifold these are domains walls, strings, monopoles and textures. Among them, 
cosmic strings and monopoles seem to be the best candidate to be observed. 
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A global monopole is a spherical heavy object formed in the phase transition of a system 
composed by a self-coupling scalar field iso-triplet ip"", whose original global 0(3) symmetry is 
spontaneously broken to U{1). The matter fields play the role of an order parameter which 
outside the monopole's core acquires a non-vanishing value. The simplified global monopole was 
introduced by Sokolov and Starobinsky [Hj. The gravitational effects of of the global monopole 
were considered in Ref. where a solution is presented which describes a global monopole at 
large radial distances. The gravitational effects produced by this object may be approximated 
by a solid angle deficit in the (3+l)-dimensional spacetime whose line element is given by 

ds^ = dt^ - dr^ - a^r^ {dO^ + sin^ ed^"^) . (1) 

Here the parameter = 1 — SvrG??^ is smaller than unity and depends on the energy scale rj 
where the global symmetry is spontaneously broken. The energy-momentum tensor associated 
with this topological object has a diagonal form and reads: Tq = = (1 — o?^lc?r'^ and 
T| = Tg = 0. The spacetime defined by the metric tensor above is not flat: the scalar curvature 
R = 2(1 — a^)/a^r^ and the solid angle is = Aira'^, smaller than the ordinary one, consequently 
there is a solid angle deficit given by 60. = 4tt{1 — a^). It is of interest to note that the effective 
metric produced in superfluid '^He — A by a monopole is described by the line element with 
the negative angle deficit, a > 1, which corresponds to the negative mass of the topological 
object 0- The quasiparticles in this model are chiral and massless fermions. 

The quantum effects due to the point-like global monopole spacetime on the matter fields 
have been considered for massless scalar jH] and fermionic fields, respectively (for the heat- 
kernels on the generalized cone see |H1 IH])- In order to develop this analysis, the scalar re- 
spectively spinor Green functions in this background have been obtained. The influence of the 
non-zero temperature on these polarization effects has been considered in for scalar and 
fermionic fields. Moreover, the calculation of quantum effects on massless scalar field in a higher 
dimensional global monopole spacetime has also been developed in ;TJ|. 

Two important ingredients responsible for vacuum polarization effects are the non-trivial 
geometry of the spacetime itself and boundary conditions obeyed by the matter fields. In this 
direction, the total Casimir energy associated with massive scalar field inside a spherical region 
in the global monopole background have been analyzed in Refs. |S1 by using the zeta func- 
tion regularization procedure. Moreover, bosonic Casimir densities induced by a single shell have 
been calculated in jl.Sj to higher dimensional global monopole spacetime by making use of the 
generalized Abel-Plana summation formula jl41ll5j . More recently, using also this formalism, a 
similar analysis for spinor fields has been developed in ^Hl- In these publications two distinct 
situations were examined: the vacuum averages inside and outside the shell. The generalized 
Abel-Plana summation formula is a powerful instrument which allows to develop the summa- 
tion over all discrete modes. By using this formula, all the physical vacuum polarization local 
quantities can be separated in two contributions: boundary dependent parts and independent 
ones. The boundary independent contributions are similar to previous results obtained in the 
literature for scalar and fermionic fields, [ElEl, using different approach. These contributions 
are divergent and consequently in order to obtain finite and well defined results, one must apply 
some regularization procedure. The boundary dependent contributions, on the other hand, are 
finite at any strictly interior or exterior points and do not contain anomalies. Consequently, 
they do not require any regularization procedure. 

Bosonic Casimir densities induced by two concentric spherical shells in the global monopole 
spacetime have been analyzed in Here we shall continue in this line of investigation, calcu- 
lating spinor densities for two concentric spherical shells. In this context, we shall admit that the 
fermionic field obeys the MIT bag boundary conditions on the spherical shells. Specifically we 
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shall calculate the vacuum expectation value of the energy-momentum tensor and the vacuum 
interaction forces between two shells. Also the total Casimir energy will be considered. 

The Casimir effect associated with fermionic fields in a Minkowski background, was con- 
sidered in a great number of papers jl8., 19^ 20 , 21j ,221 UHl US US] (for reviews and additional 
references see (22111311111231201)) considering that the fields obey the MIT bag model condition. 
As far as we know, the most of the previous studies were focused on global quantities, such as the 
total vacuum energy and stress on the surface. The density of the fermionic vacuum condensate 
for a massless spinor field inside the bag was investigated in Ref. HOI (see also |29j). As to the 
Casimir effect it is of physical interest to calculate not only the total energy but also the local 
characteristics of the vacuum, such as the energy-momentum tensor and vacuum condensate. In 
addition to describing the physical structure of the quantum field at a given point, the energy- 
momentum tensor acts as the source of gravity in the Einstein equations. It therefore plays an 
important role in modelling a self-consistent dynamics involving the gravitational field For 
a = 1 case, our calculation is a local extension of the previous contributions on the fermionic 
Casimir effect for a spherical shell. 

This paper is organized as follows: In SectionlHwe consider the vacuum expectation values of 
the energy-momentum tensor and the fermionic condensate in the region between two spheres. 
Section 131 is devoted to the investigation of the vacuum interaction forces between the spheres. 
In Section IS by using the zeta function technique, we investigate the total Casimir energy in 
the region between two spheres as a sum of the zero-point energies of elementary oscillators 
and various limiting cases are discussed. In Section [51 we present our concluding remarks. In 
Appendix El by using the generalized Abel-Plana formula we derive a summation formula for 
the series over the eigenmodes of the fermionic vacuum in the region between two spheres. 

2 Vacuum expectation values of the energy-momentum tensor 
and the fermionic condensate 

In this section we shall consider a massive spinor field propagating on a point-like global 
monopole background described by the line element . The dynamics of the field is described 
by the Dirac equation 

i7'^(9^ + r^)^-MV' = , (2) 

where 7^ are the Dirac matrices defined in such curved spacetime, and = 7i/V^7^/4 is the 
spin connection with being the standard covariant derivative operator. Here we shall adopt 
the following representation for the Dirac matrices^ 



given in terms of the curved space Pauli 2x2 matrices a 



k. 



I _ I COS 9 e"*"^ sin ^ \ 2 _ ^ f - sin 6* e"**^ cos ( 



~ arsin^ [ e^^ ' ' ^ ' 



These matrices satisfy the relation 



^i^k ^ ^ik ^ ^ p ^ 
V7 



Ikm 



^Notice that using this representation the usual anticommutation relation between the Dirac matrices, 
{7m 5 71-} = Sg^i^ is observed. 
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where 7'*^ = —g''^ are the spatial components of metric tensor and 7 is the corresponding 
determinant, e''^'" is the totally anti-symmetric symbol with e^"^^ = 1, and the latin indices 
i,k, - ■ ■ run over values 1, 2, 3. 

In this paper we are interested in the vacuum expectation value (VEV) of the energy- 
momentum tensor in the region between two spherical shells concentric with the global monopole 
on which the field obeys bag boundary conditions: 

(l + z7'^n(r'^) VL=^ = , «; = a,6, (6) 

where a and b are the radii for the spheres, a < b, n^^^ = rS^^5^ is the outward-pointing normal 
to the boundaries. Here and below we use the notations n^") = —1, n*^*) = 1. To find the VEV 
for the energy-momentum tensor we expand the field operator in terms of the complete set of 
single-particle states |V'x^\V'x 

X 

where a^^ is the annihilation operator for particles, and 6j is the creation operator for antipar- 

ticles. Now we substitute the expansion ((Tj) and the analog expansion for the operator into 
the corresponding expression of the energy-momentum tensor for the spinor fields, 

^/.i. = ^ 'l7(/.Vj,)V^ - (V(^|)7^)V; . (8) 

By making use the standard anticommutation relations for the annihilation and creation oper- 
ators, for the VEVs one finds the following mode-sum formula 

(0|V|0) = ^V{vSi-Hx),4-)(x)} , (9) 

X 

where |0) is the amplitude for the corresponding vacuum. For the problem under consideration 
the eigenfunctions have the form jlHj 

Tin 



E = y^WTM^,l = j-^,n^ = i-ir, (11) 

where a = (cr^, cj^, cr^), and n = f/r. These functions are specified by the set of quantum 
numbers x = (crkjm), where < A: < 00, j = 1/2, 3/2, . . . denotes the value of the total angular 
momentum, m = — j, ■ ■ ■ ,j determines the value for its projection, a = 0, 1 specifies two types 
of eigenfunctions with different parities corresponding to I = j — with / being the orbital 
momentum. In formula (|10)) . ^Ijim are the standard spinor spherical harmonics whose explicit 
form is given in Ref. j32j, and Zy^[x) represents the cylindrical Bessel function of the order 

yo = ^- 12 

a I 

The functions (jlUI) are orthonormalized by the condition 

d3xV7V'?+4''^ = W , ^ , V = ±, (13) 
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from which the normahzation constant A can be determined. As in mode-sum formula © the 
negative frequency modes are employed, in the discussion below we shall consider the eigen- 
functions (|lUjl with the lower sign. Note that in terms of the function Z^^ikr) the boundary 
conditions ^ take the form 

Zu„ikw) = —Z^^+n„(kw)—= ^ ^"^^ , (14) 

with w = a,b. 

In the region between two spherical shells concentric with the global monopole, the function 
Zu^{kr) is a linear combination of the Bessel, Jui^), and Neuamnn, Yy{z), functions. The 
coefficient in this linear combination is determined from the boundary condition © on the 
sphere r = a and one obtains 

4:HfcO = t]{ka, kr) ^ J,^{kr)Yj^\ka) - Y,^{kr)jj,l\ka), (15) 

where for a given function F{z) we use the notation 

F^^\z) = zF'iz) + L(-) {^Jiu, - V^i^TT^) - n^z^J F{z) , (16) 



with w = a,b, and = Mw. Now from the boundary condition on the outer sphere one finds 
that the eigenvalues for k are solutions to the equation 

C,^ (b/a, ka) Jif {ka)Yi;^ (kb) - F ^ (^a) j£) (kb) = 0. (17) 

Below we denote by 'yu^,s = ka, s = 1,2,..., the positive roots to equation (|17l) . arranged 
in ascending order, 'yu^,s < 7u^,s+i- Substituting the eigenfunctions into the normalization 
integral (|13)) with the integration over the region a ^ r ^ 6 and using the standard integrals 
for cylindrical functions (see, for instance, j33j). for the normalization coefficient of the negative 
frequency modes one finds 

TT^ki^/WTM^ - M)^ , , , , 
^ = — T^jr], ka), ka = ju„ s, 18 

where Tjy(ry, ka) is defined in Appendix 1X1 

Since the bulk and boundary geometries are spherically symmetric and static, the vacuum 
energy- momentum tensor has diagonal form, moreover {Tq) = (T^"). So in this case we can 
write: 

(0 |t;| 0> = diag(e, -p, -p^, -p^) , (19) 

with the energy density e, radial, p, and azimuthal, p±, pressures. As a consequence of the 
continuity equation V^, (O |T^| O) = 0, these functions are related by the equation 

r^+2(p-px) = , (20) 
ar 

which means that the r-dependence of the radial pressure necessarily leads to the anisotropy in 
the vacuum stresses. 

Substituting eigenfunctions (|1U() into Eq. Q, the summation over the quantum number 
m can be done by using standard summation formula for the spherical harmonics. For the 
energy-momentum tensor components one finds 



oo 

-vr 



^ (2j + l) Y.^uAri.z)fit z,g\^]{z,zr/a) 



j=l/2 (7=0,1 s=l 



(21) 
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with q = £, p, p± and we have introduced the notations 



with g]^^^^{z,y) = Ju+n„{y)Yu"'' {z) — Yy^^^{y)Jy'^'{z). The vacuum expectation values given by 
formulae H21|) are divergent and need some regularization procedure. To make them finite we 
can introduce a cutoff function ^\{z)^ z = ^ua,s with the cutoff parameter A, which decreases 
sufficiently fast with increasing z and satisfies the condition <I>a ^ !> A ^ 0. Now we apply 
to the sum over s the summation formula derived in Appendix EI As a function h{z) in this 



z-^{2u + n„) 



gi^\z,y)gii(z,y) , 



(23) 
(24) 



formula we take h{z) = f, 



the function fif)^ 



/a) ^x{z). As it has been pointed out in Appendix El 



the right of formu 



z,gi^J {z, zr/a) satisfies relation and, hence, the part of the integral on 
a over the interval (0, /x^) vanishes after removing the cutoff. As a result 

(25) 



the components of the vacuum energy-momentum tensor can be presented in the form 

= qia{r) + qab{r), 



where 



and 



j = l/2 cr=0,l-^^ -Jya 



x,g^u}{x,xr/a) 



(26) 



-1 



i=i/2 



(7=0,1 /3=± 



r "1 

+ E E / dxniftl{ax,bx)F^'fJ ax, G^f ) (ax, rx) , (27) 



with r2i^^(ax,6x) defined in Appendix 1^ In formula (|27j) we have introduced the notations 



Fif) \xMfKax,y) 



F^f) \x,G'^f\ax,y) 



F^P^f^ \x,Gl^^\ax,y) 



(Vx2 - M2 + piM)G';^'^^^{ax,y) 
(\/x2-M2 - /?iM)G(f£ (ax, y)" 

'G('^^)2(ax,y)-GK(ax,y) 



Vx2 - M2 L 



-rin 



2u + Tlo- 

y 



G("^)(ax,y)G(f] (ax,y) 



f^^C(-)(a.,.)Gffi(ax,.) 



where 



G('"±)(x,y) = /,(y)i^(-±)(x)-i^,(y)/^±)(x), u; = a,6. 



(x,y) 



(y)i^(-±)(x)+if,+„„(y)4-±)(x). 



(28) 

(29) 
(30) 



(31) 
(32) 
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In these formulae and below for a given function F(z) we use the notation F^'^^^z) defined by 
formula H93|) in Appendix The expression on the right of Eq. (|27|) can be further simplified 
by taking into account the relations 



(33) 



(34) 



where cr' = for a = 1 and a' = 1 for cr = 0. From these relations and formulae ()28 |) -()3U |) . ()92() 
it follows that 



il'^\x,y) 



-K(^^)2(x)J7W(x,y), 



X, G, 



x,y) 



{a±)2. 



kI:^^\x) 



(35) 
(36) 



^UQ (x) 

By using these formulae we see that the different parities corresponding to cr = 0, 1 give the 
same contributions to the VEVs of the energy-momentum tensor and the formula H27|) can be 
written as 



Qabir) 



-1 



oo 



1=1 i3=± 



poo p 

[ax^ rx 

a_i JM L 



where 



I 1 
a 2 



and the notation (|93|) is specified to 



(37) 



(38) 



(39) 



As it has been shown in Ref. jl6j . the term ()25|) presents the components for the vacuum 
energy-momentum tensor in the case of a single spherical shell with radius a in the region r > a. 
After the subtraction of the part coming from the global monopole geometry without boundaries, 
Qmi"!"), the components qia{r) are presented in the form 



qia 



(r) +qair), 



where the part 



Qa{r) 



-1 °° f 



dx 



Pu '{ax) ^(0) 



1=1 f3=± 



Im Ki''^\ax) 



lu 



'x,K^{rx)] , 



(40) 



(41) 



is induced by a single sphere with radius a in the region r > a. This quantity diverges on the 
boundary r = a with the leading divergence (r — a)~^ for the energy density and the azimuthal 
stress, and (r — a)~^ for the radial stress (see Ref. [IT)]). 
Note that by using the identities 



) [,,K.(.x)] = ^^i^Fif [x,Urx)] 



+ ^ n^''"^nll^}{ax,bx)Fif '> \x,G'^J"^\wx,rx) 



(42) 
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with the notation 



li''^\ax)/ll'^\bx) 



Ki'"'\ax)ir'{bx) - ir'{ax)Kr>{bx) 



rim 



(43) 



the vacuum energy-momentum tensor in the region between the spheres can also be presented 
in the form 

q{r) = qm{r) + qb{r) + qba{r), (44) 



where 



qbair) 



-1 



27r2a2 



Y^lY. I dx9!'£{axM)F[f x,G^^^\hx,rx 



(45) 



and the quantities qb{r) are the vacuum densities induced by a single shell with radius b in the 
region r < b. The formula for the latter is obtained from ()41() by the replacements a ^ b, 
I ^ K. The surface divergences in vacuum expectation values of the energy-momentum tensor 
are the same as those for a single sphere when the second sphere is absent and are investigated 
in Ref. In particular, the term qab{^) {qbai''')) is finite on the outer (inner) sphere. It follows 
from here that if we decompose the VEVs as 



qir) = qm{r) + ^ qw{r) + q^^^Xr), 



(46) 



w=a,b 



when the interference term q'^"'^^ (r) is finite for all values a ^ r ^ b. Two equivalent representa- 
tions for this term are obtained from formulae (|37|) . (|4H) . (|45|) . It can be explicitly checked that 
separate terms on the right of formula (|46|) satisfy the continuity equation (|2()|) . 

Having the components of the energy-momentum tensor we can find the corresponding 
fermionic condensate (OlV'V'lO) rnaking use of the formula for the trace of the energy-momentum 
tensor, Tn = Mipip. It is presented in two equivalent forms corresponding to w = a,b: 



mm 



{Om\lplp\Om) + W), 



1=1 i3=± 



[wx, rx) 



rx 



°° , xn^J^}(ax,bx) 
ax — — 

M Vx^ - M2 



(47) 



where the boundary-free part (Om,|V'^|Om) and the single sphere-induced part {il^ip)w are investi- 
gated in Ref. JS] • Alternatively one could obtain formula H47() applying the summation formula 
(|91jl to the corresponding mode-sum X^^'i/'x ''V'x fermionic condensate. In the limit 

a ^ with fixed values for the other parameters, the main contribution into the interference 
parts of the energy-momentum tensor and the vacuum condensate comes from / = 1 terms and 
these quantities vanish as a2/". In the limit 6 ^ oo for a massless field the interference parts 
behave like 6~2/a-2 -poi a massive field and under the condition Mb ^ 1 the main contribution 
comes from the lower limit of the integral and the VEVs are exponentially suppressed. 

The limit a <C 1 corresponds to a large solid angle deficit. In this limit and for a fixed 
r the scalar curvature is positive and large. In accordance with formula H38|) the order of 
the Bessel modified functions in the formulae for the interference parts is large. By using the 
uniform asymptotic expansions for these functions (see, for instance, j3Sj ) . we can see that the 
interference parts are suppressed by the factor exp[(2/a) ln(a/6)]. As it has been mentioned in 
Introduction, the global monopole geometry with a negative solid angle deficit corresponding to 
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a > 1 arises in superfluid ^He — A. So it is of interest to consider also the limit a ^ 1. Note 
that in the limit a — > oo one has — > 1/2 and the series over I in the expressions for the VEVs 
diverge. It follows from here that for a ^ 1 the main contribution into these series comes from 
large values I and to the leading order we can replace the summation over I by the integration: 

Y^lf{u)^a^ dy{y - 1/2) f{y), (48) 
1=1 -^1/2 

where we have used the expression for v from (|38j) . Now by taking into account formulae 
(|S7)) . (|17|) . we conclude that in the limit under consideration the boundary induced VEVs 
for the local observables tend to the finite limiting value obtained from these formulae by the 
replacement . Note that in this case the scalar curvature for the background spacetime also 
takes the finite limiting value — 

In the formulae above taking a = 1 we obtain the corresponding quantities for a spinor field 
in the Minkowski bulk. In this case u = 1 + 1/2 and the Bessel modified functions are expressed 
in terms of elementary functions. Note that the previous investigations on the spinor Casimir 
effect for a spherical boundary consider mainly global quantities, such as total vacuum energy 
or the vacuum forces acting on the boundary. For the case of a massless spinor the density of 
the fermionic condensate {ipip)h is investigated in |20| (see also |29jl. 



3 Vacuum interaction forces between the spheres 

An interesting property of the Casimir effect has always been the geometry dependence of the 
vacuum forces acting on the boundaries. In this section we investigate the interaction forces 
between the spheres as functions on sphere radii and the solid angle deficit. The vacuum force 
per unit surface on the sphere at r = if is determined by the radial stress evaluated at this 
point. The corresponding effective pressures can be presented as a sum of two terms: 

pH=p(-)+p(-)^, ^ = a,6, (49) 

where the first term on the right is the pressure for a single sphere sX r = w when the second 
sphere is absent and the term is induced by the presence of the second sphere. The 

term pf''^ directly evaluated from the corresponding bulk stress in the limit r — > is divergent 
due to the surface divergences in the subtracted VEVs and needs additional renormalization. 
This can be done, for example, by applying the generalized zeta function technique to the 
corresponding mode-sum. This calculation lies on the same line with the evaluation of the total 
Casimir energy for a single sphere in the global monopole spacetime and will be presented in 
the forthcoming paper. The second term on the right of Eq. ()49() determines the force by which 
the fermionic vacuum acts on the sphere due to the modification of the spectrum for the zero- 
point fluctuations by the presence of the second sphere and can be termed as an interaction 
force. This term is finite for all nonzero distances between the spheres and is not affected by 
the renormalization procedure. As the vacuum properties are r-dependent, there is no a priori 
reason for the interaction terms (and also for the total pressures p^"^^) to be equal for w = a 
and w = b, and the corresponding forces in general are different. For the sphere at r = w, the 
interaction term is due to the summand p^""^^ for the inner sphere and for the outer sphere. 
Substituting into these terms r = a and r = b respectively and using the relations 

G';^'"'^\wx,wx) = -1, G^^Jl\wx,wx) = n(^)5(^Va;, (50) 
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with the notation 
one finds 



1 



P(int) 2TT^a^W ^ ^ 
1=1 I3=± 



dx 



M 



xQ.wl{ax^ bx) 



M2 



wa 



(51) 



(52) 



Now by making use of the Wronskian for the Bessel modified functions it can be seen that the 
following relation takes place 



In 



1 



li''^\ax)Ki''^\hx) 



bW^ 



2v-l 



w dw L Kh^'''{ax)ir"'{bx)\ L w 
This allows us to present the vacuum interaction forces in the form 



l^(£)(ax,6x). (53) 



P{mt) 



n 



M Q 



1=1 



dx 



M 



M2 



In 



1 



ll''^\ax)Ki'+\bx) 



ax)li^^\bx) 



(54) 



where we have taken into account that F^'^~\wx) = F^'^~^^*{wx). In the next section we will 
show that these forces can also be obtained from the Casimir energy by standard differentiation 
with respect to the sphere radii. In particular, we shall see that in the limit a, 6 — > oo, with 
fixed b — a, from ()54|) the result is obtained for two parallel plates with bag boundary conditions 
in the Minkowski background. Other asymptotic expansions for the interaction forces can be 
obtained from the corresponding formulae for the interaction energy and will be discussed in 
the next section. 



4 Casimir energy 

Due to the surface divergences in the VEVs of the local physical observables, the renormalized 
total vacuum energy can not be directly obtained by the integration of the corresponding energy 
density. In this section, to consider the total vacuum energy for the configuration of two con- 
centric spheres, we follow the procedure which have been frequently used in the calculations of 
the Casimir energy for various boundary geometries and is based on the zeta function technique 
(see, for instance, j34( I35j and references therein). In the region between the spheres the total 
vacuum energy is the sum of zero-point energies of elementary oscillators: 

oo CO 

Ea<r<b = -Yl + 1) E Y.(^Ls/^' + ^')'^'- (55) 
j=l/2 o-=0,l s=l 

To regularize the divergent expression on the right of this formula we introduce the related zeta 
function 

oo 

C(n)=^"+i E(2i + 1)0M, (56) 
i=i/2 

where the parameter /i with dimension of mass is introduced by dimensional reasons, and the 
partial zeta function is defined as 

oo 

0(-)= E Y.(^Ls/a' + MY''^'. (57) 

o-=0,l s=l 
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The computation of the Casimir energy requires the analytic continuation of the zeta function 
to the value u = —1. The starting point of our consideration is the representation of the partial 
zeta function as a contour integral in the complex plane z: 



(n) 



^ / dziz^ + Mr-/^§-ln[z-^C.Av,az)], 

(7=0,1 ""^ 



(58) 



where C^^ {rj, az) is given by ((T7j) and C denotes a closed counterclockwise contour enclosing all 
zeros Ju^,s- The additional factor in the argument of the log function is introduced to cancel the 
pole at z = 0. We assume that the contour C is made of a large semicircle (with radius tending 
to infinity) centered at the origin and placed to its right, plus a straight part overlapping the 
imaginary axis and avoiding the points ziziM by small semicircles in the right half-plane. When 
the radius of the large semicircle tends to infinity the corresponding contribution into (^) 
vanishes for Re n > 1 . The integral on the right of Eq. H58|) can be presented in the form 



(n) 



^2^i ^ ^ jc« 



k(ext) 



dziz"^ + Ally's In 
oz 



j!.:\az)Ht\bz) 



Ht\az)Ji:>{bz 



(59) 



where and are the upper and lower halves of the contour C, hI^'^\z) are the Hankel 
functions and we have introduced the notations 



Cf"^ in) 



^(ext) / X 



-E 

o-=0,l 



dz(z^ + M^)'2 — ln 
c dz 



1 



27ri ^-^ ^-^ I 

(7=0,1 n=l,2 



dz (z^ + M^)~2-— In 
oz 



z'^-+^-'^mj^^)(az) 



(60) 
(61) 



As before, the additional factors in the arguments of the log functions in (|60j) and (|61|) are 
introduced in order to cancel the poles at the origin. First of all let us consider the function 
()6U() . Noting that the eigenmodes in the region inside a single spherical shell with radius b are 



zeros of the function (bz), we see that C) ^ (u) is the partial zeta function for this region. It 
can be further simplified by making use of the relation 



.(int) 



n 

(7=0,1 



M 



Jlibz) - '^Mbz)J,^i{bz) - Jl^^ibz) 



(62) 



where v is defined by formula (jSHI) with / = j + 1/2. The factor (V^^ + M'^-M) j z^ on the right 
of this relation has no zeros and singularities inside the contour C and, hence, does not give 
contribution into the zeta function. Now by making use the standard properties of the Bessel 
functions, we see that the parts of the integrals over (0, itiM) cancel and we find the following 
integral representation for the interior zeta function 



c 



(int) I 



vr 



sm ■ 



/=1 

Mb^v 



dx (x^ -M^)"^^lni X 



-2f+2 



bx 



I,{bx)ll{bx) + {l + ^+ " 



bx^ b^x^ 



Il{bx) 



.. (63) 
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By the same way it can be seen that the zeta function corresponding to the part 1)6 1|) is presented 
in the form 

oo 



c 



(ext) 



U 



2 1 . ITU 

sm — 

TT 2 



^ poo 

Vm dx{x^-M 
1=1 



ox 



Kl\ax) + 2^^^±^K,{ax)Kliax) + + ^ + Kl{ax) 



ax 



ax" 



..(64) 



Now the complete zeta function can be written in the form 

2 



C(n) 



C^'"*^ (n) + C 



(ext) 



{u) + -n 

TT 



M+1 



sm ■ 



(a; 



M 



ox 



&^\ax)Kr'{bx) 



1 



ax)/^^'^^(6x) 



(65) 



I3=±' 

where we have used relation ()33() to see that different parities corresponding to cr = 0, 1 give 
the same contribution into the interference part. The term (^u) on the right of this formula 
is the zeta function for the region outside a single sphere with radius a. This can be seen by 
various ways. First of all in the limit 6 ^ oo we expect that the interference effects should 
disappear and the zeta function should be a sum of the zeta functions for single spheres. As we 
consider the region a ^ r ^ b this sum should involve the interior zeta function for the sphere 
with radius b and the exterior zeta function for the sphere with radius a. Indeed, in the limit 
b ^ oo the last term on the right of formula (|65() vanishes and the zeta function is the sum of 
single sphere parts. Another way to see that (u) is the zeta function for the exterior region 
of a single sphere is to use the method based on the scattering theory and described in Refs. 
|3n| I35j . We could also directly write down the expression for C^^^*^ (n) by using the property 
that the quantities (both local and global) for the exterior region of a single sphere are obtained 
from those for the interior region by the replacements 1^ ^ Ki, in the corresponding formulae. 
It can be easily seen that for the special case a = 1 formulae and coincide with the 
integral representations of the zeta functions on the Minkowski bulk derived in Ref . j2l] . The 
last term on the right of formula ()65() is finite at u = —1 and the further analytic continuation 
is necessary for the single sphere zeta functions only. The corresponding procedure is standard 
and has been multiply used in the calculations of the Casimir energy (see ^] and references 
therein). For the case of a massive spinor field in the Minkowski bulk this is demonstrated in 
|24j . Note that already in this special case the calculations are lengthy. The corresponding 
results for the global monopole bulk depend on additional parameter a and will be presented 
in our forthcoming paper. Here we will be concentrated on the interference part of the vacuum 
energy. 

On the base of the representation for the zeta function, the total vacuum energy in the region 
a ^ r ^ 6 is presented in the form 



Ea^r^b = C{n)\u=-i = i^il + E^."' + AE 



(66) 



where E, 



r^b 

r>a (-^r<b) vacuum energy for the region outside (inside) a spherical shell with 



(a) 



radius a (6) and the interference term is given by the formula 



AE 



IT 



1=1 



± 



dx va 



d 



M 



M2--ln 
Ox 



li''^\ax)Ki'^\bx) 



oo 



vr 



1 



dx 



M 



M2 



:ln 



ax) I I 
lif^\ax)Kl'^\bx) 



m 



(bx) 



Ki^'+Hax] 



li'^-Hbx) 



(67) 
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where we have taken into account that F^"^ \wx) = F^'^^^*{wx). Now comparing with formula 
()54() we see that the interaction energy and the interaction forces are related by the formula 



which corresponds to the standard relation p = —dE/dV, with V being the volume. By using 
the relations 



n 



li''^\ax)Ki^^\hx) _ [xI^-i{ax)/U{ax) - ^(+)] [xK^jbx) / K^-i{hx) - 
Ki''+\ax)li'^\bx) ~ [xK,iax)/K,_i{ax) + B(-)] [xl,_,{bx) / ^bx) + Bi+) 



(69) 



and the inequalities /j^_i(x) > Ii,{x), Ky{x) > Ky-i{x), /jy_i(x)-fr^_i(x) > Ii,{x)Ky{x) for the 
Bessel modified functions, it can be seen that the integrand in Eq. (|H7|) is positive and, hence, 
the interaction part of the vacuum energy is negative. Moreover, the integrand is an increasing 
function with respect to a and a decreasing function with respect to b. In accordance with 
formula (|^5|) this leads to negative pIjjjI^ corresponding to attractive interaction forces between 
the spheres. To obtain the total Casimir energy for the configuration of two spheres, E, we need 
to add to the energy in the region between the spheres, given by Eq. (|67j) . the energies coming 
from the regions r ^ a and r ^ b. As a result one receives 

E = ^('^^ + ^C') + AE, (70) 

where E^'^^ is the Casimir energy for a single spherical shell with radius r = w in the global 
monopole spacetime. This energy can be evaluated by using the integral representations of the 
zeta functions for the regions inside and outside of a single sphere given by formulae (|63j) and 
()64() . As it has been mentioned above, the corresponding procedure is similar to that previously 
employed in Refs. [241 125j for the calculation of the fermionic Casimir energy for a spherical 
shell in the Minkowski background and will be presented in the forthcoming paper. 

Formula 1)6 7|) for the interaction part of the vacuum energy can also be obtained by the 
integration of the interference part of the energy density over the region between the spheres. 
Noting that the functions Gi^^\wx,rx) and G^^^^ {wx,rx) are modified cylindrical functions 
with respect to the argument r, we can evaluate the corresponding integrals on the base of 
standard formulae jHH] for the integrals involving the square of a modified cylindrical function. 
By making use of relations H50() and additionally 

|.G(-«(„,,,)| . + ,71) 



^G'^fl(n,x,y)\ = ^n(-)i?(/3) _ 1, (72) 

one finds 



AW / dr r\^'^') (r) = Ei%^^ + E^l^ + AE, (73) 



with AE defined by formula (jlTTj) . In (|73|) . -E'^<r<oo (-^o^r^a) vacuum energy for a single 

spherical shell with radius a (6) in the region b^r<oo{0^r^a). In deriving H73|) . we have 
used the formula for e^"'^''{r) in terms of £ab{'>^) from (|37)) at the lower limit of integration and in 
terms of eba{f) from ()45() at the upper limit. Now adding to the part H73|) the parts coming from 
single spheres, corresponding to the term g«,(r) in (|46|) . for the boundary part of the vacuum 
energy in the region a ^ r ^ 6 we obtain the formula H66|l . Hence, we have explicitly checked 
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that the interaction part of the vacuum energy evaluated from the sum of the zero-point energies 
of elementary oscillators coincides with the corresponding energy obtained by the integration of 
the local density. Note that for a scalar field, in general, this is not the case and in the discussion 
of the energy balance it is necessary to take into account the surface part of the energy located 
on the boundaries (see, for instance, Refs. ISHl for various special cases of boundary geometries 
and Ref. j37j for the general case of bulk and boundary geometries). In this case the relation 
(|68|) between the vacuum forces and the total energy is also modified by the presence of an 
additional term coming form the surface energy. 

The general formula for the interaction part of the vacuum energy is simplified in limiting 
cases. First we consider the limit when a, 6 — > oo and 6 — a is fixed. In this limit the main 
contribution comes from large values I and we can use the uniform asymptotic expansions for 
the Bessel modified functions for large values of the order (see, for instance, ||2Sj)- Replacing the 
Bessel modified functions by the leading terms of these asymptotic expansions, we introduce a 
new integration variable u = \/ — and replace the summation over I by the integration: 
Si^i /(O ~^ ^ /o°° diyf{ai'). Further, introducing polar coordinates in the {u, u) plane, after the 
integration of the angular part one obtains the following result 



AE 



1 

-— 77 / dxx\fx 



2 - M2 In 



1 I ^-2(h-a)x 

x + M 



(74) 



The expression on the right coincides with the interference part of the Casimir energy per unit 
surface for parallel plates with bag boundary conditions in the Minkowski bulk (see, for instance, 

m)- 

In the limit a ^ the main contribution to the energy comes from I = 1 term. By using the 
formulae for the Bessel modified functions for small values of the argument in the leading order 
one finds 



8M / a\i , xi 



- --t^oT-tUi dx^= (75) 



ii'\x)Kl'\x) + (x2 - f,l)Ux)K,{x) + (x2 - 



ii'^\x) + {x^ - fiimx) 

where ly = 1/a + 1/2 and for a function F{z) we have introduced notation 

F("')(x) = xF'{x) + (n^^'Vt^ + i^)Fix). (76) 
This formula is further simplified in the case of a massless spinor field: 

In the limit 6 ^ oo assuming that Mb ^ 1, the main contribution into the integral over x in 
()67() comes from the lower limit and to the leading order we have the formula 



AE 



Mb ba ^ K^-i{Ma) + K^{May 



with the exponentially suppressed interference part. In the case of a massless spinor field and 
for large values b the asymptotic behavior of the interference part in the vacuum energy is given 
by formula (f77|) . 
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Now let us consider the interaction part of the vacuum energy in the hmiting cases for the 
parameter a. The Hmit a <C 1 corresponds to strong gravitational fields. In accordance with 
formula (|38j) . in this limit the order of the Bessel modified functions is large. By using the 
uniform asymptotic expansions for these functions, we can see that the main contribution comes 
from I = 1 term and to the leading order one finds 



AE 



2V2 {a/b) 



2/a 



(79) 



In the same limit, to the leading order for the interaction forces between the spheres one has 
^(int) ~ A.E/{27raw'^). In the limit a ^ 1 corresponding to negative solid angle deficit the main 
contribution into the interaction part of the vacuum energy comes from large values /. By the 
way similar to that used in section |2 for the vacuum densities, it can be seen that to the leading 
order over a one has 



4a' 



dy 



1/2 



1/2) 



M 



M2 



In 



("+) 



{ax)Kf^\bx 



K, 



(«+) 



{ax)lf^\bx) 



(80) 



It follows from here that the interaction energy per unit surface tends to the finite value. In 
figure ^ we have plotted the interaction part of the vacuum energy as a function on the ratio 
of the sphere radii for different values of the parameter a determining the solid angle deficit. 
The full curves correspond to the massless fermionic field and the dashed curves are for the 
massive field with Mb = 1. The curves for a = 1 correspond to the bulk geometry of Minkowski 
spacetime. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



alb 



Figure 1: Interaction part of the vacuum energy as a function of a/b for various values of the 
parameter a. The full curves correspond to the massless fermionic field and the dashed curves 
are for the massive field with Mb = 1. 



5 Conclusion 

In this paper we have studied the fermionic Casimir effect for the geometry of two spherical 
shells in a idealized point-like global monopole spacetime. The VEVs of the energy-momentum 
tensor and the fermionic condensate are investigated in the region between two spheres assuming 
that on the sphere surfaces the MIT bag boundary condition is satisfied. The unrenormalized 
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expectation values are expressed as series over the zeros of a combination of the cyhndrical 
functions given by formula H17|l with the tilted notation from ()16() . To deal with this type of 
series, in Appendix ^ we derive summation formula by making use of the generalized Abel- 
Plana formula. The application of this formula to the corresponding mode-sums allowed us to 
extract from the VEVs the parts due to single spheres when the second sphere is absent and 
to present the additional part in terms of exponentially convergent integrals. In particular, the 
latter are convenient for numerical calculations. We have shown that the different parities for the 
fermionic eigenfunctions give the same contribution into the VEVs. For the points away from the 
boundaries, the boundary induced terms are unamiguously defined and the ambiguities in the 
renormalization procedure in the form of an arbitrary mass scale are contained in the boundary- 
free parts only. In particular, for a massless fermionic field the boundary induced energy- 
momentum tensor does not contain conformal anomalies and is traceless. We have discussed 
the behavior of the VEVs for the energy-momentum tensor and fermionic condensate in various 
limiting cases. In particular, for strong gravitational fields corresponding to small values of the 
parameter a, these VEVs are exponentially suppressed. For large values a corresponding to the 
negative solid angle deficit, a ^ 1, the VEVs tend to the finite limiting value. By using the 
radial vacuum stress, in section |31 we consider the vacuum forces acting on the spheres. These 
forces are presented in the form of the sum of self-action and interaction terms. Due to the 
surface divergences in the VEVs of the energy-momentum tensor, the self-action part of the 
vacuum force is divergent and needs further renormalization. This can be done by making use 
of zeta function renormalization technique. The interaction parts in the vacuum forces are finite 
for all non-zero distances between the spheres and are given by formula (|52|) . An alternative 
representation for these forces is given by Eq. (|54j) . The Casimir energy in the region between 
two spheres is considered in section [l] By using the zeta function technique we have presented 
it as the sum of the single sphere and interaction energies. The latter is given by formula (|()7j) 
and is negative. We have explicitly checked that the interaction parts in the vacuum forces and 
in the vacuum energy are related by standard formula (|68|) . In addition, we have shown that 
the interaction part of the vacuum energy evaluated from the sum of the zero-point energies of 
elementary oscillators coincides with the corresponding energy obtained by the integration of 
the local density. Note that for a scalar field, in general, this is not the case and in the discussion 
of the energy balance it is necessary to take into account the surface part of the energy located 
on the boundaries. We have investigated the interaction part of the vacuum energy in various 
limiting cases. First, we have checked that in the limit a,b ^ oo, with fixed b — a, the result for 
the geometry of parallel plates with bag boundary conditions in the Minkowski bulk is obtained. 
In the limit a ^ and for fixed values for the other parameters the interference part of the 
vacuum energy vanishes as a^^". For large values b this energy behaves like for a massless 

field and is exponentially suppressed for a massive field under the condition Mb ^ 1. For small 
values of the parameter a corresponding to strong gravitational fields the vacuum energy is 
suppressed by the factor For large values a corresponding to a negative solid angle 

deficit, the interaction part of the vacuum energy per unit surface tends to the finite value given 
by formula (|5n)) . 



Acknowledgement 

A AS was supported by PVE/ CAPES program and in part by the Armenian Ministry of Ed- 
ucation and Science Grant No. 0124. ERBM thanks Conselho Nacional de Desenvolvimento 
Cienti'fico e Tecnologico (CNPq.) and FAPESQ-PB/CNPq. (PRONEX) for partial financial 
support. 



16 



A Appendix: Summation formulae over zeros of C,^{r],z) 

As we have seen in section |21 for the massive spinor field in the region between two spheres the 
eigenmodes for the quantum number k are expressed in terms of the zeros of the function 



(81) 



where the tilded quantities are defined as ((TB|) and we will assume that r/ > 1. To evaluate the 
VEV of the energy-momentum tensor we need to sum series over these zeros. Here we derive a 
summation formula for this type of series by making use of the generalized Abel-Plana formula 
derived in |14 | I15j. In this formula as functions g{z) and f{z) we choose 



r(2b). 



Hi'''\z) Hr'{z) 
with the sum and difference 



r(2a). 



hjz) 



fiz) 



h{z) 



Hi'''\z)Hr'{z) 



(2a), 



9{z) 



-i)"7(^) 



M^\r,z) hjz) 



n = 1,2. 



(82) 



(83) 



The conditions for the generalized Abel-Plana formula written in terms of the function h{z) are 
as follows 

|/i(2;)| < ei(x)e'''^l, \z\ ^ oo, z = x + iy, (84) 

where c < 2{r] — 1) and £i{x)/x — > for x +oo. Let ^^,3 be positive zeros for the function 
Cu{r], z). To find the residues of the function g{z) at the poles z = j^^g we need the derivative 



d_ _ 4 j!)'\r^z) 

9. ^^'^'"^ 7rTf{^,z) 



(85) 



where we have introduced the notation 

Tf{rj,z) = : 

with 



Ji'^^'iz) 



w 



k'-" + {M - E){M - n^vn^'"^/w) + 



n 



2wE 



(87) 



and ka = z. By using this relation it can be seen that 



vr 



Res g{z) = -Tf{r,,^,,s)Hju,s). 
z=-iv,s 4^ 

Let h{z) be an analytic function for Rez > except possible branch points on the imaginary 
axis. By using the generalized Abel-Plana formula, similar to the derivation of summation 
formula (4.13) in JH]) it can be seen that if it satisfies condition (|84|1 . 



and the integral 



h{ze''') = -h{z) + o{z-^), z 
h{x)dx 



I 

Jo 



j(")2^x) + y) 



(a)2. 



(89) 
(90) 
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exists, then 



1- K±''(-,...)Tf(,n...) - [" 



-—Res 
2 z=o 



^0 jr\x) + Yr'{x) 

hiz)M^'\rjz) 



C,irj,z)Hi"'\z) 



l-OO 

/ dxn(£){x,vx)h{xe^-^/'), (91) 



where ^i,,m < L < ^u,m+i- Here the function ^\^J {x,rjx) is defined as 



^i^Hx ,x) = Kl^'\v^)/Kl^'\^) 



(92) 



and for a given function F{z) we use the notation 



F('"±)(z) = zF'(z) + 



Fiz) . (93) 



In section |^ formula H91() is used to derive the VEV of the energy momentum-tensor for the 
region between two spherical shells in the global monopole spacetime. As it can be seen from 
expressions (|22 [) - ()24() . the corresponding functions h{z) satisfy the relation 

hixe""'/^) = -h{xe-^'l'^) , for ^ x ^ /i^ • (94) 

By taking into account that for these values x one has F^'^+\wx/a) = F'^'"-\wx/a), we conclude 
that in this case the part of the integral on the right of Eq. ipJT|) over the interval (0, //„) vanishes. 
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